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A bstrac t 

A threc"dimcnaional, v.'iscoua computer code was 
used to calculate the mixing downutream oT a typical 
turbo fan mixer geometry ^ Experimental data were 
obtolned ueing pressure and temperature rakes at the 
lobe and nozzle exit stotions* Secondory flow vel- 
ocities wore also obtained. These data were used 
to validate the computer results. An assessment was 
also made to do tormina the relative importance of 
turbulence in the mixing plicnomcnon as compared 
with the strearawise vorticity setup by the secondary 
flows. The observations suggest that the generation 
of streamwiso vorticity appears to play a signifi^ 
cant role in determining the temperature distribu- 
tion at the nozzle exit plane. 

Nomenclature 

E energy 

h metric coefficient 

i unit vector 

Pr Prandtl number 

Q magnitude of U» equal to (u^ + + w^) 

Re Reynolds number 

U velocity vector ^ * 

u prxniary velocity 

v,w secondary velocities 

x,y ,2 three component directions 

Y specific heat ratio 

p viscosi'^y 

6 vorticity 

p density 

^ velocity potential 

stream function 
Subscripts J 

I inviscid component 

p primary component 

8 secondary component 

T turbulent 

V viscous component 

(ft Irrotational component 

^ rotation 

1,2,3 components in x,y,z direction 
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Introduction 

Turbofan forced mixer nozzles have been studied 
extenslvoly as a means of Improving fue!l efflcloncy 
and providing jet noise reduction by reducing noz- 
zle exit velocity, In the lobed mixer nozzle, the 
hot core flow passes up through the lobed section 
and the cooler fan flow passes down through the val- 
leys (fig* 1) » Very simply, these nozzles increase 
mixing by providing a much greater interface be- 
tween the hot core flow and cooler fan flow. By 
mixing the core and fan flow in this manner, a jraall 
but significant performance gain can be realized, 
The level of gain depends on the trade-offs between 
the degree of mixing of the two streams and the vis- 
cous losses incurred in the mixing process. 

Three-dimensional viscous analysis has heenv 
used by Birch, Paynter, Spalding, and Tatchell^^^ 
to Investigate the mixing process within turbofan 
engine exhaust nozzles with good results. It was 
suggested, however, that the discrepancy between 
prediction and measurement may be due in part to the 
turbulence model and to the uncertainty in the Ini"* 
tlal conditions. In addition* it was also conclu- 
ded that the calculations were sensitive to initial 
boundary layer thickness on the plug. In the Aero- 
rvv'ipulslaix 1979 Conferenced^) held /it JT^ewls Research 
Center it was further suggested that the secondary 
flows may also play an important role In the mixing 
process within engine exhaust nozzles. Calcula- 
tions were presented based on the analysis of 
Kreskovsky, Briley, and McDonald that streamwiso 
vorticity, which is o measure of the degree of sec- 
ondary flow, can be enhanced within the mixer noz- 
zle due to transverse pressure gradients, arising 
from the plug shape, d2) This suggests that three 
elements affect the turbulent mixing process, that 
is, (1) the turbulent transport process, (2) the 
streamwlse vorticity at entry, and (3) Induced sec- 
ondary flow within the mixing duct itself. 

To further investigate and illuminate the phy- 
sical processes that can be important within forced 
mixer nozzles, and to explain the large differences 
in mixing with relatively small lobe geometry chan- 
ges, an experimental and analytical effort was ini- 
tiated, This paper describes some initial results 
of this effort. 

Governing Equations 

In this study, the forced mixer flow field is 
corr/,tUed by a spatial marching method which solves 
the simplified form of the three-dimensional NS 
equations . d3-5) ^ curvilinear orthogonal coordin- 

ate system is used with coordinate directions x,y,z 
and corresponding metric coefficients hx,h 2 ,h 3 » 

Here x is defined as the streamwlse marching di- 
rection and y and z are the tv;o crossflow direc- 
tions. In order to allow the use of a marching pro- 
cedure, x-dlrection diffusion terms have been ne- 
glected. The solution method centers around the 
decomposition of the velocity into primary and sec- 
ondary flow components 
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where U repreaenta the sole component of Up, 
the Iropoacd pressure and Py(x) represents the rasan 
viscous pressure drop which is determined from th^ 
mass flux condition 



Once a solution for vorticity is obtained, a vec- 
tor potential is determined from the equation 
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Tiftc rotational components of the secondary flow arc 
thus determined from the expressions 
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To account for thermal mixing between the fan and 
core streams, an energy equation is introduced in 
the form 
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The gas law is then Introduced in the form 
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The resulting velocity field is then corrected to 
satisfy local continuity through the introduction 
of a scalar potential governed by the expression 
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The irrotational secondary velocity components can 
thus be determined from the equations 
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which relates the imposed and viscous pressure gra- 
dients to the other dependent variables. 

To determine the secondary velocity. Ug is as- 
sumed to consist of an irrotational and rotational 
component. The rotational components of Ug are 
determined from solution of a streamwise vorticity 
equation of the form 
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where 


The secondary velocities are given by the sum of the 
rotational and irrotational components 

V, - 

W, • -f 

and the vector velocity is given by 

“ ■ S“p +IjW, 

The turbulent viscosity Is determined using the two 
equations K - e turbultjnce model presented by 
Launder and Spalding. W) 

The governing equations were applied to the 
general geometry shown in figure 1. The computa- 
tional mesh starts at the lobe exit plane and the 
three-dimensional '/iscous solution is marched down- 
stream. Although the mixer is axisymmetric, the 
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ClovT threo-4imen0ionaX due to assimutJmX varia« 
tion o£ the cold fan and hot core atrcmna, 

Experlmontai Setup. 

A photograph of the teat apparatus la shown In 
figure 2* It consisted of a fixed upstream model 
section and a rotating shroud* Figure 3 shows a 
cutaway of the test apparatus* The mixer lobe sec- 
tion was interchangeable* The rotating shroud con-* 
tained instrumentation rakes for probing the flow 
field* Total temperatures were moasurcd at five 
stations in the region. The first station 

was at the lobe exit plane, the second was halfway 
to the plug end, the third was at the and of the 
plug, the fourth was midway between the plug end 
and the no?«le exit, and the fifth station was at 
the nozzle exit plane* In addition, the tempera- 
tures and pressures were measured upstream for both 
the core and fan streams. The rotating mechanism 
can be seen in figure 2 as well as the rakes at the 
lobe exit and nozzle exit stotions* Total pres- 
sures were also measured at the lobe and nozzle exit 
stations* Temperature data were token over a 54- 
dcgree segment in 3 degree increments* A typical 
temperature contour was based on data for 18 angu- 
lar positions and 14 radial positions for a total of 
252 measurements* The ratio of the shroud length 
to the inside shroud diameter was 0.71 based on the 
values at the lobe exit plane, 

The test matrix of mixer geometries Is shown in 
figure 4* penetration (ratio of lobe radius to 
shroud radius) and lobe number were varied while 
maintaining constant mixing perimeter. The spacing 
ratio, that Is, the ratio of core-to-fan lobe in- 
cluded angles*, was also changed in going from con= 
tour A to contours B and C. The two lower contours 
on the plot (A and C) were designed following cur- 
rent mixer trends* Model B was designed with a 
more gentle curvature on the upper core shape so as 
to avoid local separation. Only some of the data 
obtained with contours A and B will be presented. 

The fan and core streams were operati d with a 
total pressure ratio of unity and a total tempera- 
ture ratio of 1*35. The Mach number of the fan and 
core streams at the mixing plane (lobe exit) was 
approximately 0*45 and the by-pass ratio was 4* 

Calculated Results 
A. Temperature and Vorticlty 

The computer program described was used to 
make a sample calculation of the mixing process and 
to evaluate the program in its present state of de- 
velopment. The starting condition for the calcula- 
tion was a uniform flow that was aligned with the 
computational mesh shown in figure 1, No vorticity 
or secondary flow was present at the start of the 
computation. 

The Reynolds number, based on the hot-stream 
velocity and the shroud diameter, was 900; the Mach 
number at the inlet plane was 0.35 in both streams 
and the temperature ratio was 2. The pie-shaped 
segment shown in figure 1, covering a half- lobe and 
fan trough, was used for the computational domain. 
Symmetry was assumed for the remainder of the geom- 
etry. The computational mesh within this segment 
Was composed of 20 radial points and 10 circumfer- 
ential. A dividing streamline was used at the plug 
base since the viscous marching method will not 


compute through a separated flow region* A few of 
the computed properties have been chosen for pre- 
sentation. Figure 5 shows the total temperature 
and streamwise vorticity. The sequence of the tem- 
perature vorticity patterns reveal some Intoreatlng 
feoturea which occur during the mixing process. At 
the lobe exit plane, there is no streamwise vorti- 
city and the total temperature pattern defines the 
lobe geometry* As the solution Initially moves 
downstream, there is relatively little change in 
the total temperature* Vortices are generated on 
both sides of the radial interface between the core 
flow and the fan flow* As the liiolution progresses 
downstream, the vortices initially generated in the 
core flow dissipate and the fan vortex pair domi- 
nate* As the vortex pair intensify near the end of 
the plug* faster tiixlng appears to take place. The 
overall effect Is to enhance mixing. The vortex 
intensity and resulting enhanced mixing are prob- 
ably due to the curvature of the flow at the end of 
the plug. This suggests that the cantorbody shape 
may be important to the mixing process. Downstream 
of the plug, the streaniwise vorticity dissipates, 
The computed results show that streomwlsc vorticity, 
which was not present at the start of the calcula- 
tions, is generated in the plug region of the mix- 
ing passage* 

B, Comparison With 18 Lobe Mixer 

Calculations were performed on an 18-lobe 
forced mixer nozzle configuration using a prelim- 
inary voralon of the analysis described which did 
not have the K - e turbulence model but rather a 
simple constant viscosity model. No secondary flows 
were present in the initial conditions for the com- 
» pti^-aei-on* The computed total temperature cir^ttours* 
at the nozzle exit are presented in figure 6(a) 
whereas equivalent experimental data are shown in 
figure 6(b). In spite of the simple turbulence 
model used In this calculation, there arc remark- 
able similarities between the analysis and experi- 
ment, It is observed that the high temperature 
zone is present in both patterns, and that the cold 
fan flow extends between the hot regions, toward 
the center of the duct. The computed contours also 
begin to show a "pinching’* effect on the hot zone 
similar to that existing In the experimental data* 
The lack of better agreement in the contours is 
probably duo to the fact that the secondary flows 
present in the experimental data at the lobe exit 
hove not been accounted for In the computation. On 
the other hand, the data show some asymmetry which 
is not possible with the analysis. In addition, 
the experimental pattern shows a hot region at the 
center of the duct which may be due to recircula- 
tion at the base of the plug and/or the presence of 
a hot layer on the plug surface. The calculation 
did not have an initial hot layer under the fan 
trough or a plug boundary layer; nor did it account 
for separation and recirculation at the end of the 
plug. 

The deep radial penetration of the cold fan 
flow and the beginning of a temperature Inversion 
was not predicted by the analysis* In addition, 
the predicted **hot spots" were at a lower radial 
position than the measurements and the general tem- 
perature patterns tended to be more elongated* 

Improvements in '^he comparison of theory and 
experimental data are anticipated, based on the in- 
clusion of an Improved mixing model, and correct 
secondary flow starting conditions at the lobe exit 
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p3,dn(i instead o£ cho use of uniform iubsh*»aXlgnad 
floVt Ke»ult:s presented in the toXXwing sections 
illustrate the Importonce of understanding the 
starting or inlet flov field* 

Experi menta l ^eouJLt R 
A* Comparison of 12A jindJI.21 L Mixern 

Temperature data obtained with the 12-^lobe 
mixers > with lobe configurations A and B are shown 
in figure 7* Vour axial stations ore shown > starts 
ing at tlie lobe exit plane (fig. 7(a)) and finish** 
Ing at the uozzIq exit plane (fig, 7(d)). The lo- 
cations of the intermediate stations are given in 
the section on Experimental Apparatus. The major 
difference in contours A and B lies in the fact that 
the velocity vector at the lobe exit plonc has an 
outward radiol component and less of a tendency for 
flow separation. Contour B has a slightly larger 
penetration. These differences would not appear > 
at first glance* to be significant changes in gepra- 
etry. However* as seen in the results* the B con- 
tour hos greater penetration, the maximum tempera- 
ture at the nozzle exit is significantly lower* and 
the temperature pattern Is radically different. 

B. PiBCuasion on Observed Secondary Plows 

Two significant features of the data in figure 
7^d) can be observed* The first feature is related 
to the fact that temperature Inversion has taken 
place* whereby the hot core flow has been displaced 
radially outward, and the resulting void filled by 
the cold fan flow. In order to model this phenome- 
non it is necessary to account for the secondary 
flows which enter the mixer passage* that Is, flow, 
components normal to the grid lines at the lobe exit 
plane* The second major feature is concerned with 
the distortion that takes place in the temperature 
contours. Strong secondary flows have distorted 
the 12B pattern, resulting in a horseshoe-^shaped 
contour. These data strongly suggest that second- 
ary flow, such os the vortices shown In figure 5, 
significantly contribute to the mixing process. 

C. LDV Measurements 

Measurements of the three veloc’-v components 
have been carried out at UTRC usinr LPV tech- 
nique. The measuremeuts were made on a plane 
normal to the centerline of the mixer nozzle and 
relative to a cylindrical coordinate system. Sec- 
ondary velocity vectors for the radial and azimuthal 
directions are presented in figure 8 for a position 
just downstream of the lobe exit plane. Measure- 
ments were made at five equally-spaced angular lo- 
cations between the centerline of the fan stream 
and the centerline of the core stream. Radial in- 
flow is observed in the fan regions and radial out- 
flow in the core regions. These secondary flow data 
may explain the horseshoe-shaped temperature con- 
tours obtained with the lobe B geometry* The radial 
flow pattern Is also consistent with the presence 
of a vortex aligned with the interface between the 
fan and core flows. Some evidence of swirl is also 
observed (fig. 8) at the plug surface. In addi- 
tion, a small vortex pair could be embedded at the 
top of the core, although the data are sparse. 

Tliese secondary flow vector measurements, along 
with the temperature contours presented in the pre- 
vious section, present strong evidence for the ex- 
istence of large magnitude secondary flow in mix- 
ers. Inclusion of these secondary flows, as start- 


ing eontll cions, will be necessary in order to im- 
prove the comparison between experiment ond calcu- 
lations* Prediction of the secondary flows genera- 
ted within the lobe section is a more difficult 
problem which deserves future attention. 

Conclusions 

The computer program described in this paper 
has been used to predict the existence of large- 
Hcale secondary flow downstream of the lobe exit 
plane. This large-scale, or macro mixing, is be- 
lieved to be an important mechanism In addition to 
the conventional turbulent shear mixing by which 
the energy exchange occurs in forced mixers. It is 
now believed that properly accounting for these 
secondary flows at the lobe exit plane, through the 
use of experimental data, will lead to good agree- 
ment botwen analysis and data. Experimental evi- 
dence for the presence of secondary flow has been 
found using LDV measurements and temperature data. 
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A three-dimensional, viscous computer code was used to calculate the mixing downstream of a 
tyidcal turbofan mixer geometry. Exporlmontal data wore obtained using pressure and temper- 
ature rakes at the lobe and nozzle exit stations. Ssoondary flow velocities were ai«o obtained. 
Those data wore used to validate the’ computer results. An assessment was also made to deter- 
mine tlic relative importance of turbulence in the mixing phenomenon as compared with the 
streamwlso vortlclty setup by the secondary flows, Tlio observations suggest tijat the generation 
of stroamwise vortlclty appears to play a significant role In determining the temperature distri- 
bution at the nozzle exit plane. 
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